Practality of pomeron-exchange processes in diffractive DIS 
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By using Monte Carlo simulation of pomeron exchange model, the dependence of 
the fractal behavior of the pomeron induced system in deep inelastic lepton-nucleon 
scattering upon the diffractive kinematic variables is found rather robust and not 
sensitive to the distinct parameterization of the pomeron flux factor and structure 
function. Based on this characteristic fractal plot, a feasible experimental test 
of the phenomenological pomeron-exchange model in DESY ep collider HERA is 
proposed. 



High energy elastic and diffractive processes have long been described by 
using the phenomenology of Regge theorjtl by the i-channel exchange of mesons 
and, at high energy, by the leading vacuum singularity, i.e. the pomerom. 
Because of the ignorance of the nature of the pomeron and its reaction mech- 
anisms, there exists different kinds of approaches and parameterizations of 
pomeron dynamics in current Regge theory! And all the calculation re- 
sults concerning the collective ^spects of the diffractive processes, supJr as 
cross sectipjijpf hard diffractioiu, the distribution of large rapidity gaps, jet 
productiormm in hard diffractive processes etc. were found very sensitive to 
the distinct parameterization of the pomeron. In this respect, it is natural to 
ask the following questions: Is there a way to test and justify the pomeron ex- 
change model by using current diffractive experimental equipment (e.g. DESY 
ep collider HERA) while the criterion to the experimental measurements does 
not depend upon concrete parameterization of pomeron? If yes, what is the 
characteristic behavior of the pomeron exchange model in the expected exper- 
imental measurements? 

Having in mind that the fractal and fluctuation pattern of the multipar- 
ticle production reveals the nature of the correlations of the spatial-temporal 
evolutions in both levels of parton and hadron and is, therefore, sensitive, to 
the interact dynamics of the high-energy procesfl'El, it has been proposed!^ to 
investigate the fractal behavior of the diffractively produced system by calcu- 
lating the scaled factorial moments of the multihadronic final state. In this 
talk, we find the dependence of the fractal behavior of the pomeron induced 
system upon the diffractive kinematic variables is rather robust and not sensi- 
tive to the different parameterization of the phenomenological pomeron model 
in the deep inelastic lepton-nucleon scattering (DIS). So the characteristic plot 
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about fractal behavior of the diffractively produced system can be considered 
as a clear experimental test of the pomeron exchange model within DESY ep 
collider HERA. 

The fractal (or intermittency) behavior of the diffractively produced sys- 
tem in DIS (and also hadron-hadron collider) can be extracted by measuring 
the g-order scaled factorial moments (FMs) of the final-state hadxons excluding 
the intact proton from the incident beam, which are defined bp 

P r c , _ 1 V> (n m (n rn ~ I) . . . (n m - q + 1)) 

m—l v ' 

where, x is some phase space variable of the multihadronic final-state, e.g. 
(pseudo-)rapidity, the scale of phase space Sx = Ax/M is the bin width for a 
M-partition of the region Ax in consideration, n m is the multiplicity of diffrac- 
tively produced hadrons in the mth bin, (• • •} denotes the vertical averaging 
with the different events for a fixed scale Sx. 

The manifestation of the fractality and intermittency in high enexgy mul- 
tiparticle production refers to the anomalous scaling behavior of FMm 

F q (8x) ~ (Sx)^^" ~ M^" , as M — ► oo, Sx — > 0. (2) 

The g-order intermittency index <fi q can be connected with the anomalous frac- 
tal dimension d q of rank q of spatial-temporal evolution of high energy collisions 
aH 

d q = 4> q /(q - 1) (3) 

Pomeron exchange in Regge theory has been used to describe successfully 
the main features of the high energy elastic and diffractive process. While 
waiting for the experimental measurements of DESY ep collider HERA (and 
also hadron-hadron collider) for the above-mentioned fractal behavior of the 
diffractively producted system, let us now take a closer look at what we can 
learn from the current pomeron exchange theory. 

Pomeron factorization allows the diffractive hard scattering cross-section 
to be written as 

d A a{ep^ e+p + X) _ d 2 a(eW ^ e + X) 

dXvdtdpdCp d(3dQ 2 . W 

where x b and Q 2 are the usual deep-inelastic variables; ccp is the longitudinal 
momentum fraction of pomeron emitted by proton; and (3 = xg/xp. The first 
factor of the right-hand-side of Eq. ([|) is the pomeron flux, i.e. probability of 
emitting a pomeron from the proton; the second factor, i.e. lepton-pomeron 
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hard cross section which is assumed to be independent of negative mass-squared 
t, can be calculated in the way that 

d 2 a(eP -> e + X) f d 2 6- hard (e + parton^e + X) 

« = y^G(/3) — 2 , (5) 

if we could figure out the density G(fi') of the quarks and gluons with fraction 
/3' of the pomeron momentum. Needless to say, hard scattering cross section 
rfo"hard(e + parton — > e + X) is to be computed in perturbation theory and the 
Born-level cross section diThard is proportional to 5{/3 — (3'). But both pomeron 
flux and structure function are still main uncertainties in pomerxp-model, and 
it is even unknown whether the pomeron consist mainly of gluon£H13 or of quark 
(see last Donnachie and Landshoff's papers in Ref. 0), although measurements 
of hard diffractive scattering have been performed in both lepton-hadron and 
hadron-hadron collider. 

In addition to these theoretical uncertainties there is also a uncertainty in 
the Q 2 evolution of the parton densities of the pomeron. Numerical calculations 
using ordinary QCD evolution equation (Altarelli-Parisi or DGLAR13), and 
GLR-MCJj equation in which the inverse recombination processes of partons 
has been taken into account, turned out thau3 the Q 2 evolution of the pomeron 
structure function can be very much different depending upon whether the 
non-linear recombination term of the QCD evolution equation is included or 
not. Furthermore, depending upon the initial parton distribution at a given 
momentum scale which is unknown, the size of nonlinear term may become 
too large for the QCD evolution equation to be reliable without further, but 
also unknown, correction. By assuming both leading and subleading Regge 
trajectory, a fit according to the NLO DGLAP evolution equations to HERA 
datH of (xp, /3, Q 2 ) has favored a rather peculiar "one-hard-gluon" dis- 

tribution for the pomeron. Since what we try to pursue in this note is to find 
out whether and in which range the fractal behavior of pomeron induced sys- 
tem depends on varieties of different parameterization of pomeron, we leave 
the possible anomalous scaling behavior in the QCD evolution processes to the 
further discussioitll 

In typical kinematic region of hard diffractive processes of DESY ep collider 
HERA (say, i.e. Mx > l.lGeV, and xp < 0.1), major properties of th^diffrac- 
tive events can be well reproduced by RAPGAP generator (see, e.g. Q&EJ). In 
the following intermittency analysis of the lepton-nucleon diffractive process, 
we use RAPGAP generator to simulate the pomeron exchange processes, in 
which the virtual photon (7*) will interact directly with a parton constituent 
of the pomeron for a chosen pomeron flux and structure function. In addi- 
tion to the 0(a em ) quark-parton model diagram (7*5 — > q), the photon-gluon 
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fusion (y<7 — > qq) and QCD-Compton (7*9 — * (75) processes are generated ac- 
cording to the 0(a em a s ) matrix elements. Higher order QCD corrections are 
provided by the colour dipole model as implemented in (ARIADNE )LZl, and the 
hadronization is performed using the JETSETH. The QEELradiative processes 
are included via an interface to the program HERACLEJo. 
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Figure 1: The seoepd-order scaled factorial moments F2 getting from MC simulation of 
RAPGAP generator versus the number M of subintervals of 3-dimensional (?7, p_i_, </>) phase 
space in log-log plot, and the intermittency index <j>2 correspondingly. The different kinds of 
points denote different parameterization of the pomeron flux factor f p p(t.Xp) and structure 
function G(/3), see text for detail. 

By chosing the concrete pomeron flux factor / p p(t, xp) and the pomeron 
structure function G{0), we generate 100,000 MC events, and calculate the 
second-order factorial moments in 3-dimensional (j),p±, <f>) phase space, where 
the pseudorapidity 77, transverse momentum p± and the azimuthal angle (f> 
are defined with respect to the sphericity axis of the event. The cumulative 
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variables X translated from x = (r],p±, <ft), i.e. 



X(x) = / p(x)dx/ / p(x)dx, (6) 

were used to rule out the enhancement of EMs from a non-uniform inclusive 
spectrum p(x) of the final produced particle£3. The obtained result of second- 
order FM versus the decreasing scale of the phase space is shown in Fig.l in 
double logarithm. There exists obviously anomalous scaling behavior in the 
pomeron induced interaction, so we fit the points in Fig. 1 to Eq. (Q) with 
least square method and extract the intermittency index (f>2- In Fig. 1, we have 
also shown the Monte Carlo result of the second-order factorial moments for 
different kinds of the parameterization of the pomeron flux factor and struc- 
ture function. In Fig. 1(a), (b), and (c), we keep the pomeron structure 
function G(/3) fixed as (3G((3) = 6/3(1 — /3) but vary the pomeron flux factor 

f pF (t,x F ) as f pP (t,x F ) = %^4" 2 «-W, f pP {t,x P ) = ^J-(6.38e-^\ + 

0.424e- 3 l*l) and f pP (t,x P ) = ^[F 1 (t)} 2 x 1 p 2ap{t) respectivehj. The frac- 
tal behaviors of the pomeron induced system keep almost unchanged for dif- 
ferent flux factor. On the contrary we keep the pomeron flux factor fixed 

as f pP (t,x P ) = %^4" 2qp « in Fig.l(a), (d), (e), and (f), but vary the 
pomeron structure function as /3G(/3) — 6/3(1 — (3), (3G((3) = 6(1 — /3) 5 , 
f3GLB) = (0.18 + 5.46/3)(l - 0) and /3G(/3) = 0.077tt/3(1 - /3) respectively (see 
Rem for all these parameterisations of pomeron flux and structure function 
which are extensively used). For a given pomeron flux, the fractal behaviors 
become weaker when the pomeron become softer. In Fig. 1(d) the parton 
distribution is as soft as that in proton, the intermittency index is smallest, 
which is understandable since if the hard parton in pomeron is involved it is 
more possible to evoke jets and then the anomalous short-range correlation in 
the final-state so that the intermittency index increases, and vice versa. 

It is of the special interest to investigate the dependence of the fractal be- 
havior of the pomeron induced system upon the diffractive kinematic variables. 
We generate 500,000 events by RAPGAP Monte Carlo generator, and divided 
the whole sample into 10 subsamples according to the diffractive kinematic 
variables, e.g. xb- For each subsample, we calculate the second order scaled 
FM and the intermittency index </>2, and to see how the fractal behavior of 
the pomeron induced multihadronic final states depends upon the considered 
kinematic variable. In Fig. 2 is shown the dependence of the second order 
intermittency index <f>2 on the different diffractive kinematic variables. Since 
it is well known that the gluon density increase sharply as xb decreases in 
small-xs region, the MC result from pomeron model in Fig. 2(a) means that 
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the anomalous fractal dimension d,2 of the diffractively produced system de- 
creases with increasing gluon density, which is not inconceivable if one takes 
into account the fact (see, e.g. BO) that the effect of superposition of fractal 
systems can remarkably weaken the intermittency of whole system. 
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Figure 2: The dependence of second-order intermittency index ij>2 in RAPGAtQ Monte 
Carlo implementation upon different kinematic variables. The different shapes of points 
denote different parameterization of pomeron flux factors and the structure functions in the 
same way as that in Fig. 1. 

Obvious dependence of 4>2 on pomeron momentum fraction ip of a hadron 
and parton momentum fraction /3 of a pomeron as shown in Fig. 2(b) and 
(c) implies that, the intermittency calculated here can not be only referred 
to the hadronization processes and there should be substantial correlations 
between the fractal behaviours and pomeron dynamics. In Fig. 2(a) and (b), 
the intermittency index 02 is less than for the lower xb and xp, which 
can be imputeiLto the constraint of the momentum conservation in the high 
energy procescil. Since the Leading Proton Spectrometer (LPS) has been 
used in ZEUS detector to detect protons scattered at very small angles (say, 
< 1 mrad), which make it possible to measure precisely the square of the 
four-momentum transfer t at the proton vertex, we also showed in Fig. 2(e) 
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the t-dependence of second order intermittency index in the i-region of LPS 
detector, i.e. 0.07 < -t < 0.4 GeV 2 . To be different fr om the results of other 
kinematic variables, the fractal index for the 7*P system doesn't depend upon 
the t. 

Especially, we calculate the intermittency index shown in Fig. 2 using 
different kinds of pomeron parameterization. We denote the different shapes 
of points in Fig. 2 for the different kinds of the pomeron parameterization, 
just in the same way as that in Fig. 1. Just as mentioned above, in con- 
ventional investigation, the pomeron theory has .been used to compare with 
the data about cross, section of hard diffractioru, the rapidity distribution 
of large rapidity gapa, and jet rapidity distribution^ and jet shapeElH etc., 
where the results of the pomeron model concerning these collective nature of 
diffractive process were found very sensitive to the parameterization of the 
pomeron, and the experujaexttal data in different aspects preferred different 
kinds of parameterizatioiflflCHla. It is remarkable that the dependence of the 
intermittency index, which concerned with the inherent scaling behaviors of 
diffractive processes, upon the diffractive kinematic variables in this implemen- 
tation of the pomeron exchange model are rather robust and almost the same 
for the different parameterization of the pomeron flux and structure function! 

Having in mind that the multiplicity measurement of multihadronic pro- 
duction is available in DESY ep collider HERA and especially the multiplicity 
moments and KNO scaling behaviours of final hadrons in deep-inelastic pro- 
cesses have been already studied in recent year£3, it is feasible and urgent to 
check this characteristic fractal plot in DESY ep collider HERA. Substantial 
revision would be necessary in the manner in which we have treated diffraction 
if it should turn out that experimental measurements differ drastically from 
this characteristic plot presented here. 
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